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Abstract 

We apply the probability distribution function method to the study of chiral prop- 
erties of QCD with quarks in the exact massless limit. A relation among the chiral 
condensate, zeros of the Bessel function and eigenvalue of Dirac operator is also given. 
The chiral condensate in this limit can be measured with small number of eigenvalues 
of the massless Dirac operator and without any ambiguous mass extrapolation. Re- 
sults for SU(3) gauge theory with quenched Kogut-Susskind quarks on the 10 lattice 
are shown. 
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1 INTRODUCTION 



One can not completely understand physics of hadrons without understanding the QCD 
vacuum, which main properties are confinement and spontaneous chiral-symmetry break- 
ing, characterized by the non- vanishing chiral condensate in the massless (chiral) limit. 
Suppose the quarks are degenerate. The chiral condensate per flavor is 

(^) = ^ 7 (TrA- 1 ), (1) 

where N c = 3 is the number of colors, V is the number of lattice sites and A is the fermionic 
matrix. The trace is taken in the color, spin and position space. 

On a finite lattice, however, direct computation of the chiral condensate from a chiral 
symmetric action leads to {ipip)\ m =o = 0, even though there would be spontaneous chiral- 
symmetry breaking in the infinite volume limit. In standard lattice simulations, one has 
to add a mass term X^w/'V' to the action and calculate (ipip) at some set of non-zero 
bare fermion mass m, and then extrapolate (ipip) to the massless limit by means of some 
modeled fitting functions (e.g., linear function, polynomial or logarithmic corrections). 
Unfortunately, such a process might not be well justified, and sometimes it gives very 
different results. Here are two well known evidences. 

(a) QED in 4 dimensions. Non-compact lattice QED experiences a second order chiral 
phase transition at finite some bare coupling constant g, where the chiral condensate van- 
ishes. The critical coupling and critical index determined by am extrapolation of are 
ambiguous. Detailed discussions can be found in Ref.[l]. 

(b) The one-flavor massless Schwinger model. In the continuum, as is analytically known, 
{i>i>)cont ~ 0.16e, where e is the electric charge. Unfortunately, a more careful study [2] indi- 
cates the value of the chiral condensate computed by mass extrapolation depends strongly 
on the mass range to extrapolate, and the fitting function mentioned above. 

In Ref. [3], an alternative, named the probability distribution function (p.d.f.) method 
was proposed to investigate spontaneous chiral-symmetry breaking 1 . This method has been 
tested in the Schwinger model [3], and applied to the study of the spontaneous P and CT 
symmetry breaking [6] and theta- vacuum like systems [7] as well as the phase transition of 
SU(2) lattice gauge theory at finite density [8]. 

In this paper, we will further explore the p.d.f. method by applying it to SU(3) lattice 
gauge theory with Kogut-Susskind (KS) fermions. The rest of the paper is organized as 
follows. In Sec. 2, we elaborate the idea of the p.d.f. method and derive some relations 
between the eigenvalues of the Dirac operator and chiral condensate. In Sec. 3, we present 
the results from numerical simulations. 

1 One can find an analog in statistical physics, e.g. in the analysis of spin glasses[4, 5]. However, to 
our knowledge, it is the first time to apply such an idea to first principle quantum field theory, i.e., lattice 
gauge theory with fermions. 
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2 P.D.F. OF THE CHIRAL CONDENSATE 



Let us characterize each vacuum state by a and the chiral condensate by (ipip) a . The p.d.f. 
of the chiral condensate in the Gibbs state is defined by 



P(c) 



^2lV a 5 (c- (^)a) > 



(2) 



with w a the weight to get the vacuum state a. P(c) is the probability to get the value 
c for the chiral condensate from a randomly chosen vacuum state. If there is exact chiral 
symmetry in the ground state, P(c) = 5(c). If chiral symmetry is spontaneously broken, 
P(c) will be a more complex function. Therefore, from the shape of the function P(c) 
computed in the configurations generated by a chiral symmetric action with exact m — 0, 
one can qualitatively judge whether chiral symmetry is spontaneously broken. 

In quantum field theory with fermions, chiral-symmetry breaking is dominated by the 
properties of the fermion fields under global chiral transformation ip — * exp(iot • T'-f^ip, 
with r the generator of the chiral symmetry group. From Eq. (2), one can define the p.d.f. 
of the chiral condensate for a single gauge configuration U : 



/ 


dlf 


[dip] exp (- 


■S f )6(c-j±?X x $(xMxj) 


/ 


'dip 


[dip] exp(-S» 



Pu(c) 



where in U J2 X ~P(x)ip(x)" , summation over color, spin and position indices is implied. 
One can also define the p.d.f. for all gauge configurations 



(3) 



P(c) 



5 c- 



N C V , 



^ip(x)ip(x) 



(4) 



where the expectation value "< ... >" is computed with the integration measure associated 
with the partition function 



Z = J [dip] [dip] [dU] exp (-S) = J [dU] exp (S g + IndetA) , 



(5) 



with A the fermionic matrix. 

To study the chiral properties quantitatively, it is more convenient to employ the Fourier 
transformed p.d.f. 



Pu(q) 



dc exp (—iqc) Pu{c), 



and 



roo 

P(q) = I dc exp (— iqc) P(c). 

J —CO 

Inserting Eq. (4) into Eq. (7), one obtains 

_ roo 

P(q) — dc exp (—iqc) P(c) 

J — oo 

dip [dip] [dU] exp (—S g — Sf) 



(6) 



(7) 



1 

Z 
1 

z 



dc exp (—iqc) S I c — 



N r V 



^ip(x)ip(x) 



dip] [dip] [dU] ex.pi-S g + J2ip(x)[A 



x,y 



■ 9 r 
I 

N r V 



ip(y) 



(8) 
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Integrating out the fermion fields, Eq. (8) becomes 



P(q) = 1 J [dU] det (A - J^Lj) exp (S g ) 

= ^/ W V detA C y exp(-S 3 + lndetA), (9) 
where / is the identity matrix. Generally, a fermionic matrix A can be decomposed as 

A = mI + iT. (10) 
Denoting Xj the j—th positive eigenvalue of T, then the determinants in Eq. (9) are 

N c V/2 

det A = J] ( A ? + ™ 2 )' 
Substituting them into Eq. (9), one obtains 

i „ N c V/2 I 2 , o • M V ' 

p„) = -j m e XP( - S9+ i„de t A) n ('- ^^ 



iV c V/2 



g 2 + 2imqN c V 



(N C V) 2 (A 2 + m : 



II h .':?'.,)• (12) 



From Eq.(12), we derive relations Eq. (16) and Eq. (21) between the chiral condensate 
and the eigenmodes of the Dirac operator. 

2.1 First relation 

According to the definition Eq. (7), the Fourier transformed p.d.f. of the chiral condensate 
can also be written as 

/oo 
dc exp (—iqc) P(c) 
-oo 

= /_^cex P H, C )(5(c-^E^ > 



iq ^ - , \ 



Its derivative with respective to q at q = is 



dP{q) 



dq 



f0 -^XW- (14) 



4 



On the other hand, using Eq. (12), the derivative of P(q) with respective to q at q = is 



dq 

Comparing Eq. (14) and Eq. (15), we have 



j NcV/2 

q =o ~*\JW § AfT^)- (15) 



dP(q) .1 1 2m 



1 / NcV/2 2m \ 

In the derivation, we haven't used the specific properties of lattice fermion formulation. 
It is just the standard one in the literature. The disadvantage is that to get the chiral 
condensate in the chiral limit lim m ^ limy_ +oo ('0'0), it requires a m extrapolation, since 
< ipip > \ m=0 = on a finite lattice. As mentioned above, the result depends on the choice 
of the fitting function, in particular near the chiral phase transition. It also requires the 
calculation of all eigenvalues of the Dirac operator. When the lattice volume is large, the 
computational task is huge and not so feasible. From Eq. (16), one can also derive the so 
called Banks-Casher formula 

lim lim = lirri7r(p(A)), (17) 

which is also frequently used in the literature. Here p(A) is the eigenvalue density. This 
formula means that the eigenvalues relevant for chiral- symmetry breaking go to zero in the 
infinite volume limit. The disadvantage is that on a finite lattice, p(0) = 0, it requires 
a A — > extrapolation using some modeled fitting function. Only in the infinite volume 
limit, the number of the eigenvalues approaching zero diverges, so that p(0) ^ 0. 



2.2 Second relation 



In a theory with continuous U(l) chiral symmetry (exactly when m — 0), the vacuum is 
characterized by an angle a G [—it, 7r] and the chiral condensate can be parameterized as 
{ijjip) a = Co cos a , where Co is the amplitude of the chiral condensate corresponding to the 
spontaneously broken continuous U(l) symmetry. According to the definition Eq. (2), the 
p.d.f. of the chiral condensate is 



P(c) 



m=0 



V w a 5 (c - (ipip)a) = 77- / da 5 (c - c cos a) 

V 1 271 J-n 

VT?' CG [-00,00], 
0, c < —c , or c > c . 



The Fourier transformed p.d.f. is then 



P(q) 



m=0 



/oo 
dc exp (—iqc) P(c) 
-00 

= — / da dc 5 (c — c cos a) exp (—iqc) 

Z7T J — IT J —CO 

1 r 77 

= — / da exp (—iqc cos a) = J (qco) , 
2tt J-tt 



(18) 



(19) 
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where Jo is the zeroth order Bessel function of the first kind. From above derivation, one 
sees that the p.d.f. depends on the symmetry group of the theory. 
In the chiral limit, Eq. (12) becomes 



By computing the second derivative of Eq.(20) with respect to q, and comparing it to the 
second derivative of J (qc ), we have the following sum rule: 



which agrees with the result from chiral perturbation theory [9] and chiral random matrix 
theory [10]. The advantage is that no A extrapolation is necessary. The disadvantage is 
that all the eigenvalues of the massless Dirac operator have to be calculated as in Eq. (16). 
When V is large, the computational task is huge and not so feasible. However, very large 
lattice is required to get stable and consistent results. 

2.3 Third relation 

One can perform similar analysis for the p.d.f. of the chiral condensate for a single con- 
figuration, defined in Eq. (3). The Fourier transformed p.d.f. for this gauge configuration 
Eq. (6) is then 




(20) 




NcV/2 



(21) 




(22) 



In the chiral limit, it becomes 




(23) 



This equation is equal to zero at 



q = N c VX j (U), 



j = l,...,N c V/2. 



(24) 



Performing similar procedures when deriving Eq. (19), one has for Pjj(q) 



m=0 



Pu{q) 



Jo(qco(U)). 



(25) 



m=0 
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This equation is equal to zero at 



where z(j) is the j— th zero of Jo- in Tab. 1, the first 40 zeros of Jo are provided. 
For V » 1, Eq. (24) should agree with Eq. (26) so that 

C ° (C7) = nM(u)< ''- 1 -- 00 - (27) 

co(C/) is the amplitude of the chiral condensate for configuration U. Averaging it over gauge 
configurations with fermions, we obtain 

CU) = („(10> - $L(£). (28) 

Neither m nor A extrapolation is necessary. In the chiral-symmetry breaking phase, a 
plateau for C(j) = const, will develop, from which the chiral condensate in the chiral limit 
can be extracted. 

The relation between eigenmodes and chiral-symmetry breaking is clear: if chiral sym- 
metry is spontaneously broken, i.e., c (U) ^ 0, according to Eq. (27) and Eq. (28), Xj 
should scale as z(j)/V. In the infinite volume limit V — > oo, the eigenvalues relevant for 
chiral-symmetry breaking are those going to zero as 1/V, which is consistent with Banks- 
Casher. 

The advantage of Eq. (28) is that to extract the value of C(j) from a plateau, only 
a few smallest eigenvalues are needed for this calculation. Of course, finite size analysis 
lim m ^ hniy^oo < >= limy^oo C(j) is remained to be done, as in all approaches. 



3 RESULTS IN QCD WITH KS QUARKS 

The most interesting application of this method is QCD. Here we would like to present the 
first data for SU(3) lattice gauge theory with KS fermions, which has the action S = S g +Sf\ 



S f = 



-ir^ ReTr ^' 

x,y 



U p = U^x)U„(x + n)Ul(x + v)Ut(x), 



4 1 



ri,{x) = (-1) 



Xl+X2 + --- + Xn-1 



(29) 



where (5 = 2N c /g 2 . In the chiral limit m — 0, there exists a U(l) subgroup of the continuous 
chiral symmetry, i.e., Sf is invariant under the following transformation 



ip(x) — > exp 



\X1+X2+XZ+X4 



ip(x), ip(x) — > ip(x) exp ia(—l) 



\Xl+X 2 +X 3 +X4, 



(30) 
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All simulations are done on the V = 10 4 lattice in the quenched SU(3) case. The pure 
SU(3) gauge fields are updated using the Cabibbo-Marinari quasi-heat bath algorithm, 
followed by some over-relaxation sweeps. 100-400 independent configurations are used for 
the measurements. 

Figure 1 shows the data for (ipip) at a stronger coupling (3 = 2.5469 using Eq. (16). 
A linear function is also used to extrapolate the data of (ipip) at nonzero fermion mass 
m E [0.005,0.1] to the chiral limit. Figure 2 shows the result of C(j) for the same (3 using 
Eq. (28); there is a nice plateau for j e [15,40]. The linear extrapolation result from Eq. 
(16) is consistent with the mean value of Eq. (28) in the plateau. 

The results at a weaker coupling (3 = 5.6263 are shown in Figs. 3 and 4. One sees there 
is a little change in the slope when using Eq. (16) to calculate (ipip) and the extrapolated 
value at m — 0. In comparison, there is a wide plateau for C(j) using Eq. (28). But both 
approaches still give consistent results. 

In conclusion, we have shown how the p.d.f. method for obtaining the chiral condensate 
in the exact chiral limit works in lattice QCD with quenched KS quarks. There are several 
advantages in using relation Eq. (28): only calculations of a small set of eigenvalues of 
the massless Dirac operator are necessary; there is no need for m or A extrapolation. This 
might be an alternative efficient method for investigating the spontaneous chiral-symmetry 
breaking in lattice QCD. It would also be very interesting to see application of this approach 
to other fermion formulations, e.g., overlap fermions or domain wall fermions. 
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l 1 
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m 

Figure 1: (4>vp)/4 as a function of m for (3 = 2.5469 using Eq. (16). The dotted line is a 
linear fit of the data to the chiral limit. 
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3 




1 


2.4048 


2 


5.5200 


3 


8.6540 


4 


11.7920 


5 


14.9310 


6 


18.0710 


7 


21.2120 


8 


24.3530 


9 


27.4940 


10 


30.6346 


11 


33.7758 


12 


36.9171 


13 


40.0584 


14 


43.1998 


15 


46.3412 


16 


49.4826 


17 


52.6241 


18 


55.7655 


19 


58.9070 


20 


62.0485 


21 


65.1900 


22 


68.3315 


23 


71.4730 


24 


74.6145 


25 


77.7562 


26 


80.8976 


27 


84.0391 


28 


87.1806 


29 


90.3222 


30 


93.4637 


31 


96.6053 


32 


99.7468 


33 


102.8884 


34 


106.0299 


35 


109.1715 


36 


112.3131 


37 


115.4546 


38 


118.5962 


39 


121.7377 


40 


124.8793 



Table 1: First 40 zeros of Jq. 
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0.12 - 



(3 = 2.5469 




Figure 2: C(j)/4 as a function of j for (3 = 2.5469 using Eq. (28). The dotted line stands 
for mean value of data in the plateau. 
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Figure 3: (tpip) / A as a function of m for (3 = 5.6263 using Eq. (16). The dotted line is a 
linear fit of the data to the chiral limit. 
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(3 = 5.6263 



0.07 



0.06 



0.05 - 



C(j) 
4 



0.04 
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Figure 4: C(j)/A as a function of j for f3 = 5.6263 using Eq. (28). The dotted line stands 
for mean value of data in the plateau. 
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